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Abstract 

In this paper, we investigate a nonlinear partial differential equation, arising from 
a model of cellular proliferation. This model describes the production of blood cells in 
the bone marrow. It is represented by a partial differential equation with a retardation 
of the maturation variable and a distributed temporal delay. Our aim is to prove that 
the behaviour of primitive cells influences the global behaviour of the population. 

1 Introduction and motivation 

This paper analyses a general model of the blood production system based on a model 
proposed by Mackey and Rey in 1993 JJJ. The initial form of this model is a time-age- 
maturity structured system and it describes the dynamics of proliferative stem cells and 
precursors in the bone marrow. It consists in a population of cells which are capable of both 
proliferation and maturation. 

In this model, the period of life of each cell is divided into a resting phase and a pro- 
liferating phase (see The cells in the resting phase can not divide. They mature and, 
provided they do not die, they eventually enter the proliferating phase. In the proliferating 
phase, if they do not die by apoptosis, the cells are committed to divide and give birth, at 
the point of cytokinesis, to two daughter cells. The two daughter cells enter immediately the 
resting phase and complete the cycle. In the resting phase, a cell can remains indefinitely. 

The model in [T^ has been analysed by Mackey and Rey in 1995 [T3J[T3], Crabb et al. 
in 1996 [3 [6], Dyson, Villella-bressan and Webb in 1996 7 and Adimy and Pujo-Menjouet 
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in 2001 In these studies, the model in [T^] was simpHfied by assuming that all cells 

divide exactly at the same age. In the most general situation in a cellular population, it 
is believed that the time required for a cell to divide is not identical between cells (see the 
works of Bradford et al. [3] about experiments on mice). To our knowledge, this hypothesis 
has been given for the first time in |12| for particular cases, and only numerical studies have 
been investigated. 

In [8] and [9], Dyson et al. considered a time-age-maturity structured equation in which 
the age for a cell to divide is not identical between cells. They presented the basic theory of 
existence and uniqueness and properties of the solution operator. However, in their model, 
the division is represented by the following boundary condition 



which considers only one phase (the proliferating one), and the intermediary flux between 
the two phases is not represented in this model. 

In [TS] (1994) and (TB] (1999), Mackey and Rudnicki studied the behaviour of solutions of 
the model considered in [12], but only in the case when all cells divide at the same age. They 
obtained a first order partial differential equation with discrete time delay r and a nonlocal 
dependence in the maturity variable A(m) due to cell replication. That is. 



They gave in [16j , which is more general than [15j , a criterion for global stability in such 
equations. However, these authors considered only the special case when the term /(m, v) in 
([2]) does not depend on the maturity variable. Moreover, the nonlocal function A(m) usually 
depends on t and the condition A(to) < to, used by Mackey and Rudnicki in the example 
V{m) = rm, r > 0, and g{m) = am, < a < 1, which yields to A(to) = a^^e^'"^TO, is not 
true in the general case, for all r > 0. 

In this work, we consider a general situation, when the age at cytokinesis is distributed 
with a kernel supported on an interval [r, r] with < r < r < +oo. This yields to the 
boundary condition ([8]). 

We obtain a first order partial differential equation with a distributed time delay and a 
retardation of the maturation variable, which depends also on the time delay. In our model, 
the nonlinear part depends on the maturity variable (see Equations ([9]) and pO]) '). 

The population behaviour of the model obtained is analysed in the case 



In the interpretation of this model, N{t, 0) corresponds to the population of the primitive 
cells type, and Condition ([3]) corresponds, for a long time, to an insufficient supply of the 
primitive cells. 

We prove that Condition ([3]) guarantees, under local stability (resp. local exponential 
stability) of the partial differential equation, global stability (resp. global exponential sta- 
bility) of the system. These results are obtained under the assumption that the cells have 
enough time, in the proliferating phase (r > tq), to sufficiently increase their maturity. 




(1) 



with A : [0,1] 
me (0,1]. 




lim N{t, 0) = 0. 



(3) 
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It is also proved that, if this minimum time of rephcation r is too large, then the destruc- 
tion of the population of stem cells (cells with small maturity) affects the entire system, and 
the population extincts for a finite time (see Corollary 15. 2p . 

The paper is organised as follows. In Section[21 we present the model. We obtain a system 
of two partial differential equations with distributed time delay and a nonlocal dependence in 
the maturity variable. We give an integrated formulation of our model by using the classical 
variation of constant formula. In Section [3l we study the local stability of the trivial solution. 
In Section [U we investigate the global stability of the delay differential equation satisfied by 
the population of the primitive cells type. We use a Liapunov functional to prove the global 
stability of the trivial solution of this equation. In Section [51 we provide a criterion for global 
stability of the trivial solution in our model. This criterion is based on Condition ([3]) and 
local stability of the trivial solution. 

2 Presentation of the model 

Each cell is characterised, in the two phases, by its age and its maturity. The maturity 
describes the development of the cell. It is the concentration of what composes a cell, such 
as proteins or other elements one can measure experimentally. The maturity is supposed to 
be a continuous variable ranging from to = to m = 1 in the two phases. 

Cells enter the proliferating phase with age a — Q. We suppose that they can be lost with a 
rate 7. Every proliferating cell with maturity to divides, at mitosis, with a rate n = k(to, a). 
We assume that the mapping m ^ K{m, a) is continuous, the mapping a 1— > «;(to, a) is 
piecewise continuous, and that there exist t > such that 



This latter assumption describes the fact that a proliferating cell which has not died must 
divide before it reaches the age a = t. Therefore, the age of a proliferating cell ranges from 
a = to a = r. 

Notice that t, which represents the minimal time a cell has to spend in the proliferating 
phase before it divides, can vanish. When t > 0, cells cannot divide immediately after their 
entrance in the proliferating phase. This is probably the most realistic situation since fast 
proliferations are usually observed in pathological cases. 

All cells are supposed to age with unitary velocity, that means 



and to mature with a velocity V^(to), which is assumed to be the same in the two phases and 
satisfies V{0) = 0, V{m) > for to e (0, 1], F G C^[Q, 1] and 



K(m, a) = 



if a < T, TOG [0, 1], 



and < T < (X) such that 




for TO G [0, 1]. 



da 
'dt 




+00, for TO G (0, 1]. 
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Since vlfi) '® time required for a cell with maturity mi to reach a maturity m2, this 
last condition on V implies that a cell with maturity zero is unable to mature. For example, 
this condition is satisfied if V{m) ~ aimP , with ai > and p> \. 

m— »0 

After the division, each daughter cell enters immediately the resting phase with age a = 0. 
If the maturity of the mother cell at the cytokinesis point was to, then we suppose that the two 
daughter cells have maturity g{m) at birth. It is natural to assume that g is an increasing and 
continuous function from [0, 1] into itself. We also suppose that g e C^[0, 1) and g{m) < to 
for TO G (0, 1). Furthermore, we suppose that g~^{m) = 1 if to > g{l). 

A cell can stay its entire life in the resting phase, so the age variable ranges from = 
to a = +00. We suppose that the resting cells can die at a rate S, which corresponds to 
cellular differentiation, and that the reintroduction rate /3, from the resting phase into the 
proliferating phase, depends on the maturity and the total density of cells in the resting 
phase. 

We denote by p{t, m, a) and n{t, m, a) respectively the population density in the prolifer- 
ating and the resting phase, at time t > 0, with maturity to and age a. We define the total 
densities of proliferating and resting cells at time t and a maturity level to, respectively, by 

P{t, m) = f p{t, TO, a)da, 

and 







C + oc 

N{t,m) = I n{t,m,a)da. 







Each phase is described by a partial differential equation. Hereafter, p — p{t, to, a) and 
n = n{t, TO, a). The two conservation equations are 

t + £ + ^^Tsr^ = - + "<'• "Oi "■ 

The functions 7 and S are supposed to be positive and continuous on the interval [0, 1]. 

For the moment, f3 is supposed to be a continuous and positive function. In fact, as we 
are studying a cellular population, we could consider only non-negative solutions ; hence, to 
give a correct definition of /3, we will suppose that the mapping x t-^ /3(to, x) is constant for 
X <0, with P{m,x) — /3{m,0). For example, Mackey [TT] considered the function /3 defined 

by 

aim) 



+ b{m) 



if x > 0, 



iM' if a; < 0, 

b(m) 

for all TO G [0, 1], where a and b are two positive and continuous functions on [0, 1] and n > 1. 

System (U)-® is completed by boundary conditions modelling the flux between the resting 
phase and the proliferating phase. 

The first condition, 

p{t,m,0) = l3{m,N{t,m))N{t,m), (7) 
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represents the efflux of cells leaving the resting phase to the proliferating phase. 
The second condition, 

n{t,m,0) ^2{g-^y{m) [ K{g-\m),a)p{t,g-^{m),a)da, (8) 



describes the fact that, firstly, cells with maturity m entering the resting phase come from 
proliferating mother cells with a maturity g~^{m) and, secondly, that age of a dividing mother 
cell is between r and r and this cell divides with a rate k. 

We also consider the maturity flow tt^ : [0, 1] [0, 1], defined for s < 0, satisfying the 
ordinary differential equation 

{duAm) , , 

- ^(Tr.M), s<0, 
7ro(TO) = TO, 

for all m £ [0, 1], which represents the evolution of cells maturity to reach a maturity m at 
time from a time s < 0. tt^ satisfies 7rs(0) = and -Ksijn) G (0, 1] for s < and to G (0, 1]. 
One can remark that, for all s < and to G [0, 1], 

7rs(TO) = hr^{h{m)e^), 

where the continuous function h : [0, 1] [0, 1] is defined by 

0, for TO = 0. 



Integrating the equations (jl]) and ^ over the age, using the boundary conditions ([7]) and 
([8|) and the method of characteristics, we obtain that N{t, to) and P{t, to) satisfy, for t >t 
and TO e [0, 1], 

d d 

—N{t,m) + —(V(m)N{t,m)) = -(S{m) + p(m, N{t,m)))N{t,m) + 2(g-^y {m)x 
at am , ^ 

(9) 

k{g~^{m),a)£_{a,g~'^{m))l3{A{a,m),N{t - a, A{a,m)))N{t ~ a,A{a,m))da, 
d d 

^P{t, TO,) + —(y{m)P{t, TO,)) = -j{m)P{t, to) + /3(m, N{t, m))N{t, to) 
dt _ dm /,„N 

(10) 

fc(TO, a)^(a, m)P{T:-a(jn), N{t — a, TT-a(m)))N(t — a, 7r_Q(TO))(ia, 



where 



£,{t,m) = exp I ~ y tI^^-sCto)) + V'{TT^s(m))ds^, 



and 

fc(m, a) = K(m, a) exp ( — / K(7rs_a(m), s)(is 
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The mapping A : [0, +00) x [0, 1] R+ is defined by 

A(s,to) = ■n-s{g'^^{m)). 

A is non-increasing with respect to s and non-decreasing with respect to m. Furthermore, 
for m e (0, 1], 

A{s,m)<'m if and only if s> / . (11) 



We can notice that the solutions of Equation (|10p depend on the solutions of Equation 
([9]) whereas the converse is not true. In fact, if we know the behaviour of the solutions of 
Equation ([9]), which describes the process of replication, then we can deduce the behaviour 
of the solutions of Equation (fTO|) . Hence, we will focus our study on Equation ([9]), and 
particularly on the stability of its trivial solution. 

Let C[0, 1] be the space of continuous functions on the interval [0, 1], endowed with the 
supremum norm. We give, by using, for example, the semigroup theory, an integrated version 
of Equation ([9]). First, we recall the expression of the semigroup generated by the linear 
operator 

A : V{A) c C[0,1] ^ C[0,1], 

defined by 

V{A) = (u e C[0, 1] ; u G C\0, 1], lim V{m)u'{m) = o|, 

and 

-{5{m) + V' {m))u{'m) — V {'m)u' (m) , if m G (0, 1], 
-((5(0) + F'(0))w(0), ifm^O. 

We can easily show (one can refer to for a similar operator) that the operator A is the 
infinitesimal generator of a Co-semi-groupe {S{t))t>Q defined for ip G C[0, 1], m G [0, 1] and 
i > 0, by 

iS{t)^l;){m) = ^P{7:.t{m))K{t,m), 
where ^ 

K{t,m) = exp< — / (5(7r_s(TO)) -I- F'(7r_s(m))(is >, 



Au(m) 



for all m G [0, 1] and t > 0. 

Hence, an integrated formulation (variation of constant formula) of Equation ([9]) is given 

by 

N{t,m) = (p{T,7T_^t--){m))K{t-T,m) + G{N){t,m) - J{N){t,m), (12) 
if t >T, and 

N{t,m) ^ ip{t,m), iftG[0,r], (13) 

for an initial condition (/3 G C([0,t] x [0, 1]), the space of continuous functions on [0,t] x [0, 1], 
endowed with the supremum norm 1 1 . 1 1 , where 

G{N) {t, m) = 2j_ (^j^ C(a, -^-{t-s) {m))N{s -a,A{a, n^i^t-s) (™))) x 

/3(A(a, 7r_((_^)(m)),iV(s - a, A{a, Tr_ (m))) da^ K{t - s,m)ds, 
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where 

C(s,m) = {g-^y{m)k{g-\m),sns,g-\m)), 

and 

J{N){t,m) = K{t - s,m)/3(7r_(t_s)(rn),7V(s,7r„(4_s)(m)))iV(s,7r_(4_s)(»7i))ds, 

for t > T and m e [0, 1]. 

One can also obtain, using the same technique, an integrated form of Equation (jlOp . 

Under classical Lipschitz conditions on the nonlinear part x i— > x[3{m,x), we can eas- 
ily show, using an iterative method and the method of steps, that, for each initial data 
€ C([0,t] X [0,1]), the problem p ^ - p^ has a unique continuous solution iV^, which is 
nonnegative. 



3 Local stability 

We are interested, in this section, in the local stability of the trivial solution iV = of 
Equation © or its integrated version p^ . First, we recall a useful definition. 

Definition 3.1. Let be a subset o/C([0,r] x [0, 1]). We say that a solution N of Problem 
(T^-(W\). with an initial condition Tp, is locally stable on the set i/, /or all e > 0, there 
exists IJ, — ii{e) > such that, ifipdfl and | j^s— ^|| < /i, then, for all {t, m) G [r, +oo) x [0, 1], 
\N'('{t,m) -Nit,m)\ < e. 

This solution N is said to be locally exponentially stable on the set if there exists e > 
such that, if Lp £ VL and \\Lp — Tp\\ < e, then there exist c > and d > such that, for all 
(t,m) G [r, +oo) X [0, 1], _ 

\N'^{t,m)~N{t,m)\<ce^'^^*~'''>. (14) 

We will say locally stable (respectively locally exponentially stable j if ^ — C([0,r] x [0, 1]). 

In this section, we make the assumption that the nonlinear function x >—>■ xf3{m, x), defined 
for X G M, is Lipschitz continuous in a neighborhood of zero, for all m G [0, 1], that is, there 
exist L > and e > such that, 

\xf3{m,x) - yf3{m,y)\ < L\x - y\, (15) 

for all |a;| < e, \y\ < e and m G [0,1]. For example, if (3 is defined by then the 
mapping a; G K. i-^ x(3{m,x) is Lipschitz continuous on R and a Lipschitz constant is given 
by L = sup„g[o^i] (a(m)/&(m)). 

In the sequel, / and ^ denote the quantities 

/;= inf (S(m) + V'(m)) and C ■= f C(a)da, 
™e[o,i] 

where ({a) := sup„g[o,i] K(™'^)l C(™ja) is given in the previous section. Then we can 
prove the following invariance result, which is a first step to prove our local stability result. 
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Proposition 3.1. Suppose that L(2C + 1) < /. If e C([0,t] x [0,1]) satisfies \\ip\\ < e, 
where e is given by I115\) . then, the solution of Problem f Jji|j -i [Tij)) satisfies 

\N'P{t,m)\ < e, 

for all {t, m) e [t, +oo) x [0, 1]. 

Proof. It is clear that the assumption in this proposition imphes that / > 0. Let e > be 
given by ^ and ip e C([0,t] x [0, 1]) such that \\(p\\ < e. 
We define the sequence {Nn)n&i by 

' \ ipit,m), lite [0,r], ^ ' 

and, for n > 1, 

' \ ipit,m), if t G [0,r]. ^ ' 

Since 7r_(t_7r-)(m) < m < (?(1), then |(/3(t, 7r_(j_-)(m))| < e. Moreover, 

\K{t-T,m)\ < e~'^*~~\ 

therefore 

\No{t,m)\<ee-"^'-^^ <e. (18) 
Let T > T be given. We show by induction that for aU n G N and aU (t, m) G [r, T] x [0, 1], 

|iV„(t,m)| <e. (19) 

For n = 0, we just showed it in (flSl) . Suppose that ([19]) is true for n > 0. Then 

\Nn+i{t,m)\ < \No{t,m)\ + \G{Nn)it,m)\ + | J(7V„)(i, m)|, 

with 



and 



\G{Nr,){t,m)\ <2CeL J e"'^'""'' ds, 
\J{Nn){t,m)\<eL I e^'^'-'Us. 



So 



|A^„+i(t,m)| < e(^e-'^*-^^ + L{1 + 2C) J e-'^'-'Us 



< £fe"^(*--)+L(l + 2C)^^ 



Set 



£(2C + 1) 
a = . 
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Noting that a < 1, we obtain 

\N^+i{t,m)\ < e(e-^(*-~)(l-a) + a), 
< £((1 -a) + a). 

That is, finally, 

\Nn+i{t,m)\ < e. 

Then, is true for every n E N. 

Remarking that the sequence {Nn)ne'N converges uniformly to TV^ for {t,m) € \t,T] x [0, 1] 
(this can be proved using the method of steps), we conclude that 

\N'^{t,m)\ < e. 

for all t £ \t,T] and m £ [0, 1]. This result is true for every T > t. So, it is true for f > t 
and the proof of the proposition is ended. □ 



The result of invariancc proved in Proposition 13.11 is fundamental to obtain the local 
exponential stability of the trivial solution of Problem (fT^ . This result is presented in the 
following theorem. 

Theorem 3.1. // we suppose that L{2( + 1) < /, then the trivial solution N = of Problem 
ilS^) is locally exponentially stable. 



Proof. Let e > be given by ^ and if G C([0,t] x [0, 1]) such that ||<^|| < e, and let T > t 
be given. 

We consider the sequence defined by (fT6 |) -(fT7 | . From the proof of Proposition 13.11 it follows 
that for aU n e N and {t, m) e [t, T] x [0, 1], 

\Nnit,m)\ <e. 

Since the mapping 



1 + 2CeP^ 

is continuous on [0, /] and reaches its maximum for p = 0, there exists p G (0, /) such that 

I-P I 
L < J- — < 



1 + 2CeP-^ 2C + 1 
Using the same arguments as in the proof of Proposition l3.11 we obtain 

\No{t,m)\ < ee-'^*-~') < ee-P^*-~^ < e. 

The estimate 

|7Vi(t,m)-iVo(t,m)| < |G(7Vo)(t,m)| + |J(7Vo)(t,m)|, 
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|J(7Vo)(f,m)| < L J_\No{s,n_^t_,){m))\e-'^*-'Us, 

< eL j\-^^*-'^e-P^'-^Us, 

< eLe-^'eP^ f e'^^'P^'ds, 



and 

l-L / i~ 



\G{No){t,m)\ < 2lJ_ (^j C(a)|iVo(s-a,A(a,7r_(t_,)(m)))|(ia^e-^(*-^)ds, 



< 2eL / e" 



yields to 



Note that 



< 2QeL j e-^^^-'h-P^'-^^Us, 

< 2C,eLe-^'e^P~ j e^^-P^ds, 

\Ni{t, m) - No{t, m)| < eLe'^^eP^ 1 + 2CeP^ J e^^'P^l'ds. 



p-it^pT / ^ii-P)s^, = e-''eP^- , 

Jr I-P 

= Q(p-I)t(.P(r-t)t ^ 



= e-P(t-r) 



I-P 

1 _ e(^-P)(--*) 



Consequently 

where 



I-P 

,-p{t-' 

I-P 

\Ni{t,m) - No{t,m)\ < £i6»e-''(*-^\ 
^ 1 + 2CeP^ 



I-P • 

We show by induction that, for all n > 1 and all {t, m) € [r, T] x [0, 1], 

\N^{t,m) - iV„_i(t,m)| < e(Le)"e-^(*-^). (20) 
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Suppose that (PO)) is true for n > 1 ; then, for all (i, m) G [r, T] x [0, 1], 

\N.n+i{t,m) ~ Nn{t,m)\ < |(G(iV„) - G(iY„_i))(t, m)| + \{J{N„) - J{N.a-i)){t,m)l 
with 

|(G(iV„)-G(A^„_i))(t,m)| < 2L 



C{a)\Nn{s - a, A(a,7r_(f_5)(m))) 
— A^„_i(s - a, A(a, 7r_((_5)(m)))|da^ e~^^*~^^ds, 



and 



|(J(iY„)- J(Ar„_i))(i,m)| <LJ_ |iV„(s,^_(4_,)(m)) 

-7V„_i(s,^_(i_,)(m))|e-^(*-^)ds, 



Hence, finally, we obtain 

|Ar„+i(t,m) -7V„(t,m)| < ei^+i^e^'e-^* [l + 2Ce'= 

This means that (1^0)1 is true for n > 1. 
Then, we have, for all (i, m) e [t, T] x [0, 1], 

n 

\Nnit,m)\ < ~ N^-lit,m)\ + \No{t,m)\, 

1=1 

n 

Since < 1, we conclude that, for all t £ [r, T] and m e [0, 1], 

|iV^(i,m)|<^^e-''(*-^). 
This result is true for any T > t. So, it is true for t >t, and the proposition is proved. □ 
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4 Global stability of the population of primitive cells 

We are interested, in this section, in the behaviour of the population of immature cells. That 
is, the cell population with maturity m — Q. 

Let Nf be a solution of Problem ©-(US]) for an initial condition (f G T^{A). We know 
that 

d 

lim V{m) — N'P{t,m) = 0. 
m— >o am 

Then, the function x defined, for all t > 0, by 

:E(i)=A^'"(i,0), 
is a solution of the following delay differential equation 

-- -{Ia+ Mx{t)))x{t) +2 [ Z{a)f3o{x{t-a))x{t-a)da, (21) 




x{t) =ip{t) ■=ip{t,0), forte[0,T], (22) 

/o = 5(o) + y'(o), 

Po{x) ^ I3{0,x), foraUxeM, 

Z{a) ^ ({a, 0) = e-'"^fc(0, a), with v = 7(0) + ^'(0). 
It is clear that Iq > 0. 

Since tp G C([0,r] x [0, 1]), then ip £ C[0,r], where C[0,r] is the space of all continuous 
functions on [0,t]. Under classical Lipschitz conditions on the function /3o, one can prove the 
existence and uniqueness of solutions of Equation ([2T|l (see Hale and Lunel [10 ). 

As in [11] , it is natural to add the assumption that the function (3q is decreasing for a; > 
and satisfies 

lim Po{x)=0. 

X — * + oo 

A first property of the solutions of Problem ([2T|l - ((22|) is given in the following proposition. 

Proposition 4.1. Let ip e C[0,t]. // '0 non-negative on [0,r], then the solution of 
Equation \21]) for the initial condition ip, denoted by x^ , is non-negative on [r, +00). 

Proof. Suppose that the assumption of the proposition is false. First, suppose that there 
exist to £ [t,t + r) and e > such that to -\- s < r + r, x^{to) — and x^{t) < for 
t G {to, to+e). Then, for aU t G [to, to + e), 

-{Io+f3o{xHt)))x^{t)>0, 
and _ _ 

Z{a)j3o{x'^'{t- a))x'^{t- a)da= / Z{a)l3o{-ip{t - a))i}{t - a)da > Q, 



because, for all t G [t,t + r] and all a G [t,t], t — a £ [0,r]. Hence, for t £ [to, to + e), 
{x'^y{t) > 0, and we obtain a contradiction. So > on [r, t + t]. 

Using the same arguments and the method of steps, we obtain the result for any t >t. □ 
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As we are studying a cellular population, it is more interesting to consider non-negative 
initial conditions. Then we introduce the set 

C+ = {Vj e C[0,r] ; > 0, for i G [0,t]}. 

The global stability of the population of primitive cells is given by the following result. 

Theorem 4.1. Assume that Iq > (2z — l)/3o(0) where z = Z{a)da. Then the trivial 
solution X = of Equation 1121]) is globally stable on the set , that is, for all G , 

lim x'^{t)=0. 

t — ^-t-oc 

Proof. We prove the global stability of the trivial solution of pTjl by using a Liapunov 
functional. 

We define the functions A and A by A(a;) = xj3(^{x) and K{x) — / X{s)ds, for all a; > 0. Let 

Jo 

H : C+ ^ M be the mapping defined, for all 'tp G C+, by 

H{i^) = K{iP{t)) + I Z{a)( f \''{ij{9))d9)da. 



We show that H is a. Liapunov functional relative to Equation (pT|) . We have 
Hi^P) = limjiHixt)~HW), 

= kr)XWr)) + 1^ Z{a){X^{i:{T)) - X\^/j{t - a)))da. 

Since 



iP(t) = -{lo + /3o(^(t)))^(t) + 2 J Z(a)A(^(T - a))da, 

then 

H{^) = -[/o+/5o(V(t))]/3o(^(t))V'2(^) 

Z(a)[A2(V.(T)) + 2A(^(t))A(V'(t - a)) - X\ij{T - a))]da, 



lo + l3oiHrmoWr))r{r) + 2 J Z{a)X\Mr))da 
Z(a) [A(^(t)) - A(V'(t - a))] ^da. 



Hence, 

Hii:) < ~[Io ~ (2z - l)/3o(^(T))]/3o(^(r))^2(-)^ 
Let u be the function defined, for x > 0, by 

u(x) — r{x)/3Q{x)x'^ , 
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where, for all x > 0, 

r(.T) =/o - (2z- l)/3o(x). 

r satisfies r(0) > and liuix^aorix) = /o > 0- Since /So is decreasing on r is a 
monotonous function. Therefore, r is a positive function on R"*". 

Hence u is non-negative on [0, +oo) and u{x) = if and only if x = 0. Consequently, every 
solution of Problem ((2T|) - ((22|) with -0 G C"*" tends to zero as t tends to +oo. □ 

Combining this result with the result of local stability shown in Section [31 we can obtain 
global stability for a class of solutions of Equation ([9]) or its integrated version ([12]). 

5 Global stability 

In this section, we give a general criterion for global stability of Equation p^. We first recall 
a definition. 

Definition 5.1. Let il be a subset of C{[0,t] x [0,1]). A solution N of Problem ifIg|j - iiTg)) . 

with initial condition Tp, is said to be globally stable on the set Q if, for all ip £ Q and all 
me [0,1], _ 

lim \N'^{t,m) -N{t,m)\ = 0. 

This solution N is said to be globally exponentially stable on the set f2 if for all ip Cz ^ there 
exist c > and d > such that, for all (t, m) G [t, +oo) x [0, 1], 

\N'^{t,m)~N{t,m)\ <ce-'^^'-^\ 

We will say globally stable (respectively globally exponentially stable j ifil — C([0, t] x [0, 1]). 

From the result of local stability obtained in Section [31 we are able to establish a first 
result about global stability for the trivial solution of Equation ([72]) . This result is presented 
in the following proposition. 

Proposition 5.1. Suppose that the mapping x ^ x[3{m,x) is Lipschitz continuous on 
[0, +cxd) for all m G [0, 1], with a Lipschitz constant L, such that 

i(2C + l)</. (23) 

Then the trivial solution N = of Problem is globally exponentially stable. 

Proof. We first show, as in the proof of Proposition [2111 that for all (p G C([0, r] x [0, 1]) and 
{t,m) G [0,+oo) X [0,1], 

\Nf{t,m)\ < \\ip\\. 

Then, we have (using the same arguments as in the proof of Theorem 13. 1|) that, 

lim N^{t,m) = exponentially, 

t — >oo 

for aU G C([0, r] x [0, 1]) and m G [0, 1]. □ 
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The inequality 

L(2C + !)</, 

is satisfied \f I or v are large or if L is small. Recall that / and v represent the mortality 
rates in the resting and the proliferating phase and that L is a bound of the quantity of 
reintroduced cells. Then, the result of Proposition 1 5 . 1 1 represents a normal behaviour of the 
cellular population. 

Corollary 5.1. // the function a; G K i-^ /3(m, x) is differentiate and decreasing on (0, +oo) 
for all m G [0, 1], then, under the assumption 

(2C + 1) sup /3(m,0) </, 

me[0,l] 

the trivial solution N = of Problem is globally exponentially stable. 

Proof. If the function x ^ x) is differentiable and decreasing on (0, +oo) for all m G 
[0, 1], we obtain that 

d 

— {x(3{m,x)) < (3{m,x) < /3(m, 0), 

for X £ (0, +oo) and m G [0, 1]. Then, the fmiction x t-^ x[3{m, x) is Lipschitz continuous on 
M and the proposition 1 5 . II concludes the proof. □ 



For example, if /3 is defined by ([6]) and if, in addition, 

(2C+1) sup </, 
me [0,1] \ Km)J 

then the trivial solution = of Problem (|12p is globally exponentially stable. 

Our next objective is to prove the global stability of solutions under more general condi- 
tions, less restrictive than the condition ((23)) given in Proposition l5.ll 

We need the following assumption : there exists C > such that 

for all m G (0, 1]. 

The biological interpretation of this condition is easy. That is, the time ^fg) 
required for a daughter cell to reach the maturity g~^{m) of its mother is bounded. For 
example, if V{m) ~ aim, with ai > 0, then ([24)1 is equivalent to g'{0) > 0. 

m— »0 

Then we can prove the following fundamental theorem. 

Theorem 5.1. Suppose that the condition ^24^ is satisfied and that the function x i— !■ 
xP{m,x) is Lipschitz continuous for all m G [0, 1]. Set 



To = sup 

m>0 



ds \ 
V{s))- 
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Suppose that Ni{t,m) and N2(t,m) are two solutions of Problem (T^-(W\) for the initial 
conditions ipi and </?2 respectively. If 

L > To, 

and if there exists b G (0,1] such that 

ipi{t,m,) = (p2{t,m), 
for m G [0, b] and t G [0, t], then there exists t > T such that 

Ni{t,m) = N2{t,m), 

for m G [0, 1] and t >t. 

Proof. Since r > tq, then, from (fTT|) . A(a, m) < m for all a G [z, t]. 

We first show, using the Gronwall's lemma and a method of steps, that Ni{t,m) = 
N2it, m) for m G [0, b] and t > 0. 

Let m G [0, b] and t > t (if t G [0, r], the result is obvious). Then 

7r_((-r) ("T") < m < b. 

So 

'Pi{T,Tr-(t-r){m)) = (^2(T,7r_(t_^)(m)). 

Let t e [t, T + t] . Then, 

A(a, 7r_(j_^)(m)) < 7r_(t_^)(m) < m < b, 

and consequently 

G(7Vi)(i,m) = G(7V2)(t,m). 
As the mapping x i— > x/3(m, x) is Lipschitz continuous and A' is bounded, we can write 

\Ni{t,m)-N2{t,m)\ = \J{Ni){t,m) ~ J{N2){t,m)\ 

where c is a nonnegative constant. Then, using the Gronwall's lemma, we obtain that, for 
all i G [t, r + r] and m G [0,6], 

Ni{t,m) = N2{t,m). 

Using the method of steps, we obtain this result for t >t. 
We define the mapping a : [0, 1] — + [0, 1] by 

a(m) ~ A(t, to). 

Notice that a is an increasing function such that a{m) < to for all to G (0, 1] and a(.g(l)) = 
h^^{e^-). From pTjl and since r > tq, /i^^(e^-) < .g(l). Hence, a is invertible and we denote 
by A its inverse function on [0, ft-^^(e^-)]. We set A(to) = g(l) for m G [h~^{e^-), g{l)]. Since 
a is non- negative and increasing on the interval [0, 1], A is increasing on [0, h^^(e^-)]. 
We first suppose that b G (0, /i^^(e^-)]. 
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We consider the increasing sequences (6„)„gN and (<„)„gN defined by 

bn+i = A(6„), for n G N, 

and 

t„+i = ^„ + T + lnf ^^^"\^^ y fornGN, 



io = 0. 



We show by induction that, if m G [0, 6„] and t > tn, then Ni{t, m) = -/V2(^, "^). 

Suppose the equahty is satisfied for n > 0. Let m G [0, bn+i] and t > tn+i- Then t > t„ + t. 

The solutions Ni, i — 1,2, of Problem can be rewritten as 



7V,(t,m) = 7V,(i„ + r,7r_(i_t„_7)(m))A'(t-i„-T,m)- J(7V,)(t,m) 
2 ^* 



( y C(7r-(t-s)(TO),a)iVi(s - a, A(a,7r_(t_s)(TO)))x 
/3(A(a,7r_(f_^)(TO)),7Vi(s - a, A(a, 7r_(t_s) (m))))da^ A'(t - s,m)ds 
for t > t„ + r. Since 

t — i„ — T > In 

we have 



h{bn) 



h{bn+i) 
Hence 

h{bn+i) h{bn+i) 
because h is an increasing function. Consequently, 

Using the induction hypothesis, we obtain that 

Nl{tn +T,Tr_(^t-t„-T){m)) = N2{tn + T,7r„(i_f„_7)(m)). 

Finally, since s — a > tn and 

A(a,7r_(f_5)(m)) < A(a,m) < A(t, m) < A(t, 5„+i) = 5„, 

we obtain that 

\Ni{t,m) - N2it,m) \ = |J(iVi)(t,m) - J{N2){t,m)\. 

Using the Gronwall's lemma, we obtain the equality for n + 1. 

Noting that there exists M G N* such that 5a/ < (?(1) = &a/+i, we obtain, for all 
m G [0,5a/+i] = [0,1] and t > tM+i, that Ni{t,m) = N2{t,m). 

We now suppose that b G {h~^{e~-), g{l)]. Then the above result is still valid with M = 0. 
Thus, 

Ni{t,m) = N2{t,m) for m G [0,5(1)] and t > tM+i, M G N. (25) 
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We finally suppose that h £ (5(1), !]• 

Let m G [3(1), 1] and t > tu+i + 
overliner be given. For s G [Im+i +T,t] and a G [l,t], we have s — a > tM+i- Thus, if 
x(-(t-s),m) <g(l), 

A{a,x{-{t~s),m)) < A{a,g{l)) < h-\e-^) < g{l) 

and, from (|25p. 

Ni{s - a, A(a, x(-(i - s), m))) = A^2(s - a, A(a, x(-(t - s), m))), 
and if x(— (t — s), m) > g{l) then, from the definition of C, 

C(x(-(i-s),m),a) = 0. 
Moreover, since /i(to) < 1 for to G [0, 1], we obtain, for to G [5(1), 1] and t > tM+i +t — 

Xi-it - tM+i - T),m) < h-\h{m)h{g{l))) < 5(1). 

We deduce that 

Ni(tM+i +T,x{-{t - tm+i -T),m)) = N2{tM+i +T,x{-{t - tM+i ~T),m)). 

Using once again the Gronwall's Lemma, we get 

Ni{t,m) = N2(t,m) for m G 1] and t > tAf+i + t - ln(/i(g(l))), M G N. 

We set t ~ <A/+i + r — ln(/i(g(l))) and the theorem is proved. □ 

In the sequel, we suppose that Condition \2A\ is satisfied, that the function x ^ x[3{m, x) 
is Lipschitz continuous for all to G [0, 1] and that t > tq. 

Corollary 5.2. Let ip G C([0, r] x [0, 1]). // there exists b G (0, 1] such that 

Lp{t, to) = 0, 

for TO G [0, b] and t G [0, t], then there exists t > t such that 

Nf{t,m) = 0, 

/or TO G [0, 1] and t >t. 

The result in CoroUarv lS. 21 describes the influence of stem cells (which are small maturity 
cells) on the total population. If the time of replication is large enough (r > tq), then a 
destruction of the population of stem cells affects all the system and the total population of 
cells dies out after a finite time t. Biologically, this case corresponds to the aplastic anemia, 
a disease which destroys or injures stem cells, preventing the division. 

Let if G C([0, r] x [0, 1]) and b G (0, 1]. We define the function ipb : [0,t] x [0, 1] ^ M by 

, . / ip{t,m), ifTOG[0, 6], , . 

^(t,5), ifTOG[M(l)]. ^''^ 

It is clear that ipb G C([0,r] x [0, 1]) for all b G (0, 1]. 

To prove our main result on the global stability, we need two preliminary lemmas. 



18 



M. Adimy, F. Crauste 



Global Stability of a PDE with distributed delay 



Lemma 5.1. Let ip S C([0,t] x [0, 1]). Then, for allT > t and e > 0, there exists b G (0, 1] 
such that 

\N'^''it,m) - N'^it,0)\ < e, 

for te[T,T + T] and m G [0, 1] . 

Proof. Let (p G C([0,t] x [0,1]), T > t and e > be given. Using the continuity of the 
solutions of Equation (fT^ . we deduce that there exists 6 = S{e) > such that, if the function 
V' G C([0,r] X [0,1]) satisfies 

m,m)-ipit,0)\<S, 

for {t,m) G [0,r] x [0, 1], then 

\N'^{t,m) - m{t,0)\ < e, 
for (t, m) G [T, T + t] x [0, 1]. Moreover, (p is continuous, so there exists b G (0, 1] such that 

\ip{t, m) — (p{t, 0)1 < S, 
for all {t,m) G [0,t] x [0,6]. Hence, we obtain that 

\ipbit,m)-ip{t,0)\ <5, 

for {t,m) G [0,r] x [0, 1], and so, 

\Nf"'{t,m) - Nf{t,0)\ < e, 
for te[T,T + T] and to G [0, 1]. □ 
Lemma 5.2. Let (p G C([0, r] x [0, 1]). Then, for all b G (0, 1], there exists t > such that 

N'^''{t,m) = N'f{t,m), 

for t >t and m G [0, 1] . 

Proof. From the definition of ipb and Theorem 15. 11 the lemma is obvious. □ 

Now, we are able to prove that the global stability of the solutions of depends on 
their local stability and the global stability of the population of immature cells. 

Theorem 5.2. Suppose that the trivial solution N = of Equation il^) is locally stable 
(resp. locally exponentially stable). Then this trivial solution is globally stable (resp. globally 
exponentially stable) on the set 

= ((p G C([0,t] X [0,1]) ; lim N'f(t,0) = o\. 

Proof. = is a locally stable solution of Equation (fT2|) . Then, for all e > 0, there exists 
S > such that, if G C([0, r] x [0, 1]) satisfies 

IMI<<^, 
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then 

\N'f(t,m)\ < e, 

for t > T and m G [0, 1]. 

Let e > be fixed and ip E fl. Then hmt^oo N'^{t, 0) — 0. So there exists T > t such that, 

\N'^{t,0)\ < d/2, 

for t >T. Furthermore, from Lemma [O] there exists b € (0, 1] such that 

\N'^''{t,m) - N'^{t,0)\ < S/2, 
for {t,m) e [T,T + t] X [0,1]. Hence, 

]N'^''{t,m)\ < 5, 

for {t,m) G [T, T + r] x [0,1]. Since Equation is autonomous, A'"'^'' becomes an initial 
condition for on the interval [T, T + t]. So, we obtain that 

|7V^''(i,m)| < e, 

for t > T + r and m £ [0,1]. 

From Lemma [5?2l there exists Z > such that 

N'^''{t,m) ^ N'^{t,m), 

for t >t and m G [0, 1], so 

\N'^{t,m)\ < e, 
for t > max{T + T,t} and m G [0, 1]. Consequently, 

lim N'^{t,m) = 0, 

t — ^oo 

for all (fi G ft and m G [0, 1]. 

If A'' = is locally exponentially stable, the conclusion is obtained in the same way. □ 

The theorem 15.21 stresses the influence of immature cells in the behaviour of the total 
population. That is, if the population of immature cells is stable, it implies that the total 
population is also stable. 

Remark that the assumption 

L(2C + 1) </, 

given in Proposition 15.11 implies, from Theorem 13.11 the local exponential stability of the 
trivial solution of Equation (fT^ . Furthermore, since /3o(0) < L, z < ( and I < la, we have 

(2z - l)/3o(0) < lo. 

Hence, it follows, from Theorem 14. 11 that 

c+ c n. 

Then, Proposition 15. II becomes a simple consequence of Theorem 15.21 

We conclude with a general result of global stability, concerning the entire system of blood 
production. 

If the trivial solution iV = of Equation (fT2)) is globally stable on a set fl, then the trivial 
solution {N,P) = (0,0) of the integrated version of Problem ([9|)- p0)) is globally stable on 
the set n X C([0,t] x [0,1]). 
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